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Abstract 

We study a two-dimensional cellular automaton (C A) , called Diffusion 
Rule (DR), which exhibits diffusion- like dynamics of propagating patterns. 
In computational experiments we discover a wide range of mobile and sta- 
tionary localizations (gliders, oscillators, glider guns, puffer trains, etc), 
analyze spatio-temporal dynamics of collisions between localizations, and 
discuss possible applications in unconventional computing. 
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1 Introduction 



In our previous studies on minimal cellular automaton (CA) models of reaction- 
diffusion chemical system we constructed [4J a binary-cell-state eight-cell neigh- 
borhood 2D CA model of a quasi-chemical system with one substrate, state 0, 
and one reagent, state 1. In that model chemical reactions were represented 
by semi-totalistic transition rules. Every cell switches from state to state 1 
depending whether sum of neighbors in state 1 belongs to some specified inter- 
val or not. A cell remains in state 1 if sum of neighbors in state 1 belongs to 
another specified interval. 

From 1296 cell-state transition rules, we selected a set of rules with com- 
plex behavior 4 . Amongst the complex rules, namely in G-class of morpho- 
logical classification [4J, we located so-called Diffusion Rule. CA governed by 
this rule often exhibits slowly non-uniformly growing patterns, resembling dif- 
fusive patterns in chemical systems with non-trivial coefficients of diffusion, or 
reaction-dependent diffusion coefficients, so the name of the rule. 

The rule simulates sub-excitable ^36. medium-like mode of perturbation 
propagation — cell in state takes state 1 if there are exactly two neighbors 
in state 1, otherwise the cell remains in state 0, and, conditional inhibition — 
cell in state 1 remains in state 1 if there are exactly seven neighbors in state 1, 
otherwise the cell switches to state 1. 

In present paper we are trying to answer the following questions. Is there 
a reaction-diffusion binary-state CA that express complex dynamic? Can we 
demonstrate that CA exhibits non-stationary growth of reaction-diffusion pat- 
terns? Do stationary or mobile generators of localizations, glider guns, exist in 
binary-state reaction-diffusion CA? Can the reaction-diffusion CA simulate an 
effective procedure and therefore be universal? 

CA with space-time dynamics similar to that in spatially extended chem- 
ical systems are studied from early days of CA theory and applications [55] , 
however most rules discovered so far lack minimality (some of the rules employ 
dozens of cell-states). Methods of selecting the rules also widely vary depending 
on theoretical frameworks, e.g. probabilistic spaces [201 138] and genetic algo- 
rithms [mini]. Therefore, we envisage a strong need for a systematic analysis 
of propagating patterns like those observed in the Diffusion Rule. The propa- 
gating patterns are of upmost importance in modern computer science because 
such patterns play a vital role in developing novel and emerging computing 
paradigms and architectures, particularly collision-based computing fTl [3| l23]. 

We must mention that various authors have already obtained pioneering re- 
sults in the studied rule. Magnier et al. discovered three primary gliders [55] . 
David Eppstein found four gliders known, already incorporated in our frame- 
work, and four new gliders which were novel for us (Fig.|4](q) (t) and (u))Q The 
ghder traveling along diagonals of the lattice (Fig.|4](v)) was firstly recorded by 
Amling in 2002 (see Eppstein's web site) . Finally, a glider gun and three puffer 
trains were discovered by Wotowiczj^ 



^ See Eppstein's findings at http : //f ano ■ Ics .uci ■ edu/ca/r ules/b2s7/ 



• ^http : //hw .mirwoj ■ opus ■ chelm. pl/ca/rules/lif 6.2 ■ glf 
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Diffusion Rule CA is just one of many complex GJ^ exhibiting mobile lo- 
calizations. Other famous examples include semi-totalistic rules as the Game of 
Life [milQ], Brain's-brain and Critters rules [34], High Life [9\, Life 1133 |22], 
Life Without Death [TS], and Life variant i?35/S'236j^ Other variant is with 
Larger than Life (TB] and the Beehive and Spiral rules hexagonal CA [Ml 151 15]. 
more recent candidates were proposed by George Maydwell with Hexagonal Life 
and Hexagonal Long Life rulesj^ Amongst 3D binary state CA supporting glid- 
ers Life 4555 and Life 5766 by Carter Bays [T] |8] are most widely known. In ID 
there are Rule 110 [TS] EH IBUl [26] and Rule 54 [HI [211 [25] CA, which support 
an impressive range of mobile localizations. 

Our paper is structured as follows. In Sect. 2 we introduce basic concepts 
of CA model under investigation. Section 3 introduces results of statistical 
analysis of the Diffusion Rule using mean field theory. In Sect. 4 we present basic 
structures discovered in the Diffusion Rule CA. Section 5 compiles a catalogs 
of non-trivial interactions between mobile localizations, which could be used to 
designing basic elements of collision-based computers. In Sect. 6 we highlight 
our achievements in analysis of the Diffusion Rule and prospects for future 
studies. 

2 Basic notations 

We study family of 2D binary-state cellular automaton (CA) defined by tuple 
(Z^, E, u, /), where Z is the set of integers, every cell x G Z"^ has eight neighbors, 
orthogonal and diagonal (i.e. classical Moore's neighborhood) u{x) = {?/ G Z : 
X ^ y and |a; — y| < 1}, E = {0, 1} is the set of states, and / is a local transition 
function defined as follows: 



where a* = \{y E u{x) : y*" — 1}|, and 6*1, <^i, ^2 are some fixed parameters 
such that < 6*1 < 02 < 8 and < (5i < (52 < 8. 

We can write the rule as R{5i520i02) or like 56*1 . . . 92/S5i . . .62 more tra- 
ditional code. Also, the rule can be interpreted as a simple discrete model of 
a quasi-chemical system with substrate '0' and reagent '1', where [6*1, 6*2] is an 
interval of reaction, or association between substrate and reagent, and [5i,(52] 
is an interval of dissociation. The family of rules includes Conway's Game Life, 
when intervals [5i,(52] and [611, ^2] are interpreted as intervals of survival and 
birth, respectively. 

^http : //uacom p.uwe ■ ac .uk/genaro/otherRules .html| 
^http : / /www. ics .uci . edu/~eppstein/ca/b35s236/ 

^Several interesting eandidates in ticxagonal representation are proposed by Andrew 
Wuensche using DDLAB at http: // www.ddlab ■ com/| and by Maydwell using SARCASim at 
[http: //www. collidoscope . com/ca/ 




1, if (a:* = and cr* 
0, otherwise 




a* e [61,62]) 



(1) 
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In our previous paper |1] we morphologically classified all 1296 rules, and 
studied how changes in parameters R{SiS20i92) of cell-state transition rule in- 
fluence space-time dynamics. For example, we discovered [4] a small subset of 
rules Life 2c22j^2 < c < 8, which could be interpreted as quasi-chemical precip- 
itating systems. For parameter set [611, 02] = [22] [UIJi the system is transformed 
into 2+-medium, CA model of excitable system in sub-excitable mode. 

We have found a cluster of semi-totalistic rules supporting structures of the 
Diffusion Rule. They are B2/S2 . . .8 called Life dc2^ where d and c take 
values between 2 and 8, and d < c. Therefore, we found that the rule B2/S7 
or i?(7722) exhibits most reach dynamics of localized patterns amongst all the 
rules studied by us. Rules of the local transition are simple: 

1. A cell in state will take state 1 if it has exactly two neighbors in state 
1, otherwise cell remains in state 0. 

2. A cell in state 1 remains in state 1 if it has exactly seven neighbors in 
state 1, otherwise cell takes state 0. 

3 Mean field approximation 

Mean field theory is a proved technique for discovering statistical properties 
of CA without analyzing evolution spaces of individual rules PHI [Ml [H] • The 
method assumes that elements of the set of states S are independent, uncorre- 
lated between each other in the rule's evolution space. Therefore we can study 
probabilities of states in neighborhood in terms of probability of a single state 
(the state in which the neighborhood evolves), thus probability of a neighbor- 
hood is the product of the probabilities of each cell in the neighborhood. Using 
this approach we can construct mean field polynomial for a semi-totalistic evo- 
lution rule [22] as follow: 

pt^i = E ( " : M p^'ff-''-' + E ( " : M pur^ (2) 

where n represents the number of cells in neighborhood, v indicates how often 
state 1 occurs in Moore's neighborhood, n — v shows how often state occurs in 
the neighborhood, pt is a probability of cell being in state 1 , qi is a probability 
of cell being in state 0. 

On the basis of outcomes of computational experiments we can suggest in- 
tervals of extreme densities of initial random configurations which leads to the 
emergence of localizations in the Diffusion Rule. In the lower limit best densities 
d are 0.004 < d < 0.015 and in the upper limit they are 0.992 < d < 0.997 for 
the first 15-20 steps of evolution (Fig. fTl. CA starting its evolution in random 



e 


http : //uncomp.uwe . ac.uk/KenEiro/Diffusion_Rule 


/lif e_2c22.html 




http : //uncomp . uwe . ac . uk/genaro/Lif e_dc22 . html 
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configuration with lower density of 1-states exhibit stationary or mobile self- 
localizations (like gliders or oscillators) at the beginning of the evolution, how- 
ever in many cases collisions between mobile localizations leads to catastrophes, 
when 1-state patterns spread all over the lattice. Random initial configurations 
with higher (upper limit) density of 1-states produce either vanishing reactions 
between localizations or symmetrical growing patterns emerged as unions of two 
or more gliders. 




Figure 1: Three random initial densities for the Diffusion Rule: (a) 0.004, (b) 
0.013 and (c) 0.995 respectively, on lattices of 200 x 200 to 18 generations. 

Thus, the mean field polynomial for the Diffusion Rule is following: 

Pt+i = 8pUt + 28pWt (3) 

The fixed point is 0.236 that represents configurations with large density of 
1-states emerging from any random initial condition (we should note that the 
fixed point for Conway's Game of Life is 0.37); this represents global density 
of 1-states necessary for evolution dynamics to stabilize. Also, we can see an 
unstable fixed point 0.05 (Fig. |2|, that implies the existence of regions with 
unpredictable behavior or complex dynamic [28| . 

Thus, p = is a super-stable point, although it is quite close to the unstable 
point. The super-stable is important, as it means a quiescent substrate, i.e., 
the state where live all the structures. Since mean field theory is just an ini- 
tial approximation, it ought to be worthwhile to gather up some Monte Carlo 
approximations for a few generations just to see if the first estimate maintains 
itself, more or less. 

We must also mention that the probability to find 'interesting' behavior is 
very low, about 0.05. Perhaps, this may be the reason why the Diffusion Rule 
was not studied before - when observing evolution from random configuration 
one more likely (with probability 0.3) to encounter a catastrophe (e.g. when 
placing three cells anywhere in Moore's neighborhood, just not in one line) then 
stable mobile localization. 
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Figure 2: Diagram of mean field polynomial for the Diffusion Rule. 

4 The Diffusion Rule Universe 

In present section we uncover a range of basic structures, stationary and mobile 
localizations, generators of localizations and polymer- like structures formed of 
the mobile localizations. 

4.1 Mobile self-localizations 

In computational experiments with the Diffusion Rule CA we discovered 26 mo- 
bile self-localizations — gliders or particles — traveling orthogonally or diago- 
nally in the lattice. Properties of the gliders, including volume, speed, direction 
of motion are listed in Table [H 



(a) (b) (c) (d) 



Figure 3: Configurations of minimal, or primary, gliders in the Diffusion Rule: 
(a) gl ghder, (b) g2 glider, (c) g3 glider and (d) g4 glider. 

Configurations of minimal gliders and compound gliders are shown in Figs.|3] 
and|4] respectively. From Table [l] we can see that 96% of gliders move orthog- 
onally, and that g23, g24, g25 and g26 are the largest gliders in the family of 
mobile localizations. First column in the table gives names of gliders. Second 
column in Table [l] represents glider's volume calculated as number of cells oc- 
cupied by the glider. Third and fourth columns are translation and period. 
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Fifth column shows the speed=c/period, where c is the maximum speed. Sixth 
column is the weight that represents the number of cells with state 1 within 
glider's volume. The last column indicates whether or not glider moves along 
columns and rows, or diagonals. 
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Table 1: Properties of gliders in the Diffusion Rule CA. 



There are two types of gliders — primary and compound |38l 126] : a primary 
glider can not be decomposed into smaller mobile localizations, a compound 
glider is made of at least two primary gliders. 

4.2 Oscillators 

The Diffusion Rule CA exhibits five types of stationary localizations known 
as oscillators. The two most common flip-flops and three blinker patterns are 
shown in Fig. [5j 

Flip-flop configurations are ol and o2 oscillators both structures flipping at 
45°. Blinkers of period four are o3, o4 and o5 oscillators. Table [2] shows basic 
characteristics of each oscillator. 

Some of the oscillators and their assembles act as eaters, i.e. stationary 
patterns that annihilate gliders colliding to them (see examples of collisions in 
Sect. 5.3). 
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(a) (b) (c) (d) (e) 



(f) (g) (h) 



(i) (j) (k) 



(1) 



(m) (n) 



(o) 



(P) 



(q) 



(t) 



(u) 



(v) 



Figure 4: Twenty two compound gliders in the Diffusion Rule CA: (a) g5, (b) 
g6, (c) g7, (d) g8, (e) g9, (f) glO, (g) gll, (h) gl2, (i) gl3, (j) gl4, (k) gl5, (1) 
gl6, (m) gl7, (n) gl8, (o) gl9, (p) g20, (q) g21, (r) g22, (s) g23, (t) g24, (u) 
g25 and (v) g26 gliders, respectively. 
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(a) (b) (c) (d) 



(e) 



Figure 5: Five oscillators in the Diffusion Rule CA: (a) ol and (b) o2 are flip- 
flops, (c) o3, (d) o4 and (e) o5 are blinkers, respectively. 



oscillator volume period weight 



01 4 2 2 

02 9 2 3 

03 10 4 4 

04 16 4 4 

05 30 4 8 



Table 2: Properties of oscillators in the Diffusion Rule CA. 

4.3 Avalanches 

Avalanches are novel structures, have not described before in related studies. 
They are assembles of adjacent gliders that cause explosive growth of rhomboid 
shaped patterns with deterministic edges and quasi-chaotic interior. Avalanches 
can be constructed from various compositions of adjacent gliders. 

Figure [6] shows an avalanche produced in composition of two gl gliders, ad- 
jacent at 90°; this avalanche pattern grows diagonally inside the third quadrant. 
The minimal volume of an avalanche is 4 x 4 with eight cells in state 1. 

One can use even number of gliders to construct symmetrical avalanches, 
two examples are shown in Fig. [7] 

4.4 Puffer trains 

A puffer train is a mobile localization which generate (leaves traces) of stationary 
localizations along its motion path. There are 16 known types of stable puffer 
trains (which produce oscillators) in the Diffusion Rule CA. Basic properties of 
puffer trains are shown in Table [3] 

A particular case of puffer train is shown in Fig. [9] This puffer bears frag- 
ments of g4 glider. All configurations of the puffer are displayed in Fig. |9] and 
there we can see that configurations shown in Fig. |9] (c) (d) (e) and (f ) can be 
interpreted as spaceships. 

Moreover, the Diffusion Rule CA exhibits dozens of non-stable puffer trains. 
In the Fig. [lO]we see five non-stable puffer trains, which produce asymmetrically 
growing, or quasi-chaotic, patterns. All discovered non-stable puffer trains have 
speed 1/c and period four. 
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(a) (b) (c) 

Figure 6: Two gl gliders produce an avalanche pattern: (a) initial condition, 
(b) configuration after 15 steps of evolution and (c) configuration after 200 steps 
of evolution. 





(a) (b) 

Figure 7: Construction of the syninietrical avalanches: (a) symmetrical 
growth initiated by two g4 gliders, configuration at 237th step of evolution, 
(b) avalanche pattern with non-trivial internal symmetries produced by assem- 
bly of two g2 glider and two g3 gliders, configuration at 237th step of evolution. 
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Table 3: Characteristics of puffer trains in the Diffusion Rule CA. 



(a) 



(b) 
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Figure 8: Fifteen puffer trains observed in evolution of the Diffusion Rule CA: 
(a) pl, (b) p2, (c) p3, (d) p4, (e) p5, (f) p6, (g) p7, (h) p8, (i) p9, (j) plO, (k) 
pll, (1) pl2, (m) pl3, (n) pl4 and (o) pl5 puffer trains, respectively. 



4.5 Mobile glider guns 

Glider guns are localized patterns that periodically lose their stability and give 
birth to traveling mobile localizations, gliders. In computational experiments 
we discovered twelve types of mobile glider guns in the Diffusion Rule CA. 

Basic parameters of the twelve glider guns are shown in Table [4] and gun's 
configurations in Fig. |11[ The most remarkable feature is that all primary gliders 
can be produced by glider guns. Some glider guns can generate two types of 



gliders at once, thus gun7 (Fig. 11") generates g2 and g3 gliders at the same 
time, however both gliders travel coupled in pairs. 

The generator gunl2 can produce one or two g4 gliders, see examples in 
Fig. |12[ This gun is also extendable, the extension is determined by the num- 
ber of ol oscillators (which should be more then two). Positions of oscillators 
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(a) (b) (c) 



(d) (e) (f) 

Figure 9: Special pl6 puffer train similar to spaceships. 



(a) 



(b) 



(c) 



(d) 



(e) 



Figure 10: Some non-stable puffer trains. 



12 



gun 


produce 


volume 


translation 


period 


speed 


weight 


move 


gunl 


gl 


30 


4 


4 


c/1 


7 


orthogonal 


gun2 


g4 


50 


4 


4 


c/1 


9 


orthogonal 


gun3 


g4 


50 


4 


4 


c/1 


9 


orthogonal 


gun4 


g4 


60 


4 


4 


c/1 


9 


orthogonal 


gunS 


g4 


60 


4 


4 


c/1 


9 


orthogonal 


gun6 


g4 


72 


4 


4 


c/1 


15 


orthogonal 


gun7 


g2 A g3 


72 


4 


4 


c/1 


16 


orthogonal 


gunS 


2xg4 


80 


4 


4 


c/1 


14 


orthogonal 


gun9 


2xg4 


90 


4 


4 


c/1 


14 


orthogonal 


gunlO 


gl A g4 


143 


4 


4 


c/1 


15 


orthogonal 


gunll 


2xgl 


154 


4 


4 


c/1 


24 


orthogonal 


gunl2 


(2xg4) V g4 


176 


4 


4 


c/1 


19 


orthogonal 



Table 4: Characteristics of glider guns in the Diffusion Rule CA. Second column 
of the table shows what type of glider each glider gun produces. 



determine number of glider streams generated by gun 12. 



(a) (b) (c) (d) (e) (f) (g) (h) (i) (j) 



(k) 



Figure 11: Configurations of glider guns in the Diffusion Rule CA. 

So far we did not find glider guns which produce diagonally-moving streams 
of gliders, neither guns generating compound gliders or stationary guns. 

4.6 Glider gun and puffer train 

There is at least one special mobile structure that combines in itself properties 



of both glider gun and puffer train. Figure 13 shows glider gun producing gl 
glider and ol oscillator each 4th step of CA evolution. This puffer-gun moves 
orthogonally, has a volume of 70 cells and weights 12 units. 
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(a) (b) (c) 

Figure 12: Extendable glider gun in the Diffusion Rule CA. 
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t=l 
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Figure 13: Configurations of puffer-gun in the Diffusion Rule CA. 



4.7 Avalanche gun 

Avalanche gun is another remarkable example of mobile generators (Fig. 14 1. 
The mobile gun produces an avalanche every 4th step of CA evolution. 



IN) 



t=i 



t=2 



t=3 



t=4 



<3 :■ 



(a) 



(b) 



Figure 14: Configuration of an avalanche gun in the Diffusion Rule CA. (a) 
production of two gl gUders to 90°, (b) development of quasi-chaotic reaction 
that destroys next avalanche conserving the avalanche gun. 

However, life-time of the gun producing each avalanche is short: when 
avalanche produced by the gun it grows and then destroys the next avalanche 
produced. 



5 Collisions between localized patterns 

The Diffusion Rule CA combines high-degree unpredictability with enormously 
rich dynamics of collisions between mobile and stationary objects. This section 
studies the outcomes of the collision reactions. 
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5.1 Forming diffusing patterns by collisions 



Gliders colliding in the Diffusion Rule CA can produce an explosively grow- 
ing diffusion-like pattern, the diffusive patterns do usually have non-stationary 
boundaries and they exhibit quasi-chaotic internal dynamics. 



i t 



(a) 



(b) 



Figure 15: Examples of diffusion-like patterns produced in collision between 
gliders in the Diffusion Rule CA. (a) pattern produced at 260th step of evolution 
after coUision between g4, g2 gliders and g3 ghder, (b) pattern produced at 260th 
of evolution after collision between two g4 gliders and two gl gliders. 



Two most distinct examples are shown in Fig. 15 In first example (Fig. 15 
(a)), g4 glider collides with g3 and g2 gliders, diffusion pattern produced is 



'lead' by three gliders and puffer train. In second example (Fig. 15 (b)), two 



g4 gliders collides with two gl gliders. The reaction produces multiple gliders, 
puffer trains, oscillators, and even vertical glider guns during their collision 
dynamics. 



5.2 Reactions between propagating patterns 
5.2.1 Soliton-like reaction 

In certain initial conditions gliders collide similarly to solitons, namely they 
restore their structure and velocity vector after collisions. In Fig. [16] we can 
see snapshots of the head-on collision dynamics between two g4 gliders (begin 
at even distance before collision). When the gliders collide they temporarily 
lose their stability, produce varieties of transient structures. After few steps of 
evolution the gliders are restored and transient structures are annihilated. 

At the moment, there is only a reaction soliton-like perhaps some others 
examples exist with more bigger gliders but initially have not more examples. 
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Figure 16: Soliton-like reaction between two g4 gliders. 



5.2.2 Eater reaction 

Eaters are stationary localizations which destroy gliders colliding into them. 
The most simple eater is built with ol oscillator. This eater destroys gl gliders 
but is shifted four cells along the glider's initial direction of motion. Both gl and 
g2 gliders are destroyed when they collide with the eater built of o3 oscillator. 
In Fig. 17 'a) you can see an eater made of two o3 oscillators destroying gl and 
g4 gliders. A 'universal' eater, which destroys all types of primary gliders is 
illustrated in Fig. [iT^b) . 



(a) (b) 

Figure 17: Examples of eater configurations in the Diffusion Rule CA. 

Using eaters one can control glider streams emitted by glider guns, thus in 
Fig. 18 we can see how the eater eliminates gl gliders produced by gunl. There 
is at least one mobile eater of gliders, this is g24 glider that eats gl gliders. 



t=0 


t=l 


t=2 


t=3 


t=4 


t=5 


t=6 



Figure 18: Eater destroys gliders emitted by glider gun. 



5.2.3 Delay reaction 

The basic delay reaction can be implemented when g3 glider collides with g4 
glider and changed to g4 glider in the result of collision, and the original g4 
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(a) (b) 

Figure 19: Delay reaction in the Diffusion Rule CA: (a) initial position of gliders 
before collision, (b) final result of delay operation, the glider traveling west 
delayed and translated southward. 

glider delayed for two time steps but also translated southward one cell per 



collision (Fig. 19) 



5.2.4 Multiplication and reduction reactions 



(a) (b) (c) 

Figure 20: Multiplication of gliders: (a) multiplication 1x1^4 reaction of 
gliders in binary collision, (b) reduction 4x4^4 reaction in glider collision, 
(c) conservative collision reaction. 



A typical reaction of glider multiplication is shown in Fig. 20 (a): two g4 
gliders are involved in head-on collision, with odd distance between glider heads 
before collision, four new g4 gliders are produced in result of the collision. Re- 
duction is implemented as multiplication of gliders, where gliders in the multi- 
plied columns are in proximity of each other. As shown in Fig. |20] (b), when we 
collide two rows of g4 gliders, four gliders in each row (eight gliders in total are 
involved in the collision), then just four new gliders are produced. Adjusting 
distance between gliders in colliding columns of gliders we can achieve almost 



any (but odd) result of multiplication (Fig. 20 (c)). 

Using multiplication reactions we can also construct arbitrary packages of 
gliders. For example, to construct a stream of packages, six gliders per package, 
we collide stream of four-glider g4 packages traveling west with a pair of g4 



gliders traveling east (Fig 21 (a)). Sequentially, all four-glider packages are 



transformed to six-glider packages (Fig. 21 (b)), the operation is halted by pair 
g3 gliders. 
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(a) 



(b) 



Figure 21: Constructing packages of g4 gliders by multiplication reaction, 
(a) initial configuration, (b) final configuration, recorded after 186th steps of 
evolution. 



5.2.5 Reflection reaction 

We discovered nine types of reflection-type collisions between stationary and 
mobile self-localizations. Let us discuss some examples shown in Fig. [22] When 
two gl gliders collide with each other (head-on collision, even distance, with 
slight shift between gliders along south-north axis) two g4 gliders are generated; 
these g4 gliders move in the direction perpendicular to original trajectories of 
colliding gliders (Fig. [22] (a)). Glider colliding with ol oscillator is reflected at 
the angle 90°, as shown for collision of g2, g4 and g3 gliders with ol oscillator 



(Fig. 22 (b) (c) (d)). 



t=0 t=7 



t^O t^6 t=6 



(a) (b) (c) 

1=0 1=9 



t=0 t=6 



t=0 t=5 



(d) (e) (f) 

Figure 22: Collisions leading to reflections: (a) two gl gliders collide with each 
other, (b) g4 glider collides with ol oscillator, (c) g2 glider collides with ol 
oscillator, (d) g3 glider collides with ol oscillator, (e) gl glider collides with o5 
oscillator, (f) g3 glider collides with o3 oscillator 

Gliders can also be derived in addition of reflections, when colliding with 
stationary localizations. Thus, when a gl glider collides with a o5 oscillator, two 
gl gliders (going in opposite directions to each other) are generated and follow 



trajectories perpendicular to original trajectory of collided gl gliders (Fig. 22 
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(e))j^ Similarly, when g3 glider collides with o3 oscillator, two g4 gliders are 
produced (Fig. [|2] (f)). 



5.2.6 Annihilation reaction 



(a) (b) (c) (d) (e) (f) 



(g) (h) (i) (j) (k) 

Figure 23: Initial positions of colliding gliders and oscillators leading to annihi- 
lation reaction, (a)-(f) binary collisions, (g)-(k) multiple collisions. 

Significant amount of collisions between localizations in the Diffusion Rule 
CA leads to annihilation of colliding patterns. Few examples of initial config- 



urations of colliding objects (leading to annihilation) are shown in Fig. 23 for 



binary collisions between gliders (Fig.p3^(a)-(f)) and multiple collisions between 



ghders and oscillators (Fig. 23 (g)-(k)) 



5.3 Computation in the Diffusion Rule 

Basic operations necessary to implement a functionally complete set of logical 



gates can be derived from collision dynamics presented in Fig. 22 Following 
paradigms of collision-based computing [2] we encode logical Truth by pres- 
ence of a glider or an oscillator, while absence of mobile or stationary objects 
corresponds to logical Falsity. 



Namely, Fig. 22 (b) (c) (d) demonstrate that stationary localizations, os- 
cillators, can play a role of mirrors thus defiecting gliders from their original 
trajectory. The mirrors can be used to route signals. Signals can be deleted, 
erased by placing eaters along trajectories of gliders, representing the signals. 
Signals can be also delayed in collisions with other gliders or oscillators. 



Collisions used to construct fanout gate are shown in Fig. 22 (e) (f), a 
glider collides to stationary localization, and two new gliders are produced in 
result of the collision. 



Dynamics displayed in Fig. 22 (a) shows a typical collision gate, where inputs 
X and y are represented by trajectories of ghders traveling East and West, re- 
spectively. While trajectories of new gliders, traveling South and North, encode 



^This reaction can synchronize multiple collisions as you can see in our example FANOUT. rie 
available from http : //uncomp .uwe . ac . uk/genaro/Dif f usion_Rule/dif f usionLif e .html 
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value of X and y. Constant Truth is made up of ceaseless stream of gliders, 
generated by glider guns. 



5.3.1 Constructing a memory device 



A memory in the Diffusion Rule CA can be constructed using basic interactions 
between gl gliders and ol oscillator, as illustrated in Fig. [24] A bit of informa- 
tion is represented in the memory unit by shaded 2x2 cells square in Fig. 24 
The bit can be read by sending a gl glider to the memory unit (top configura- 
tions in Fig. 24 glider travels East). When gl glider collides with oscillator ol 
oscillator, both glider and oscillator are annihilated (and one new ol oscillator 
is formed at some distance from the memory unit). Then the bit is erased as a 
result of the reading operation. To write down the bit one can send another gl 



glider (in Fig. 24 this 'writing' glider travels West) toward now empty memory 
unit and associated ol oscillator. Both gl glider and associated oscillator are 
destroyed, however ol oscillator is restored in memory unit (shaded region in 
Fig. 24 1, i.e., we write a bit again. 
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Figure 24: Constructing a memory device in the Diffusion Rule CA. Snapshots 
of the configurations. Time increases from right to left and from top to bot- 
tom. The domain where the bit of information (ol oscillator) is written to is 
represented by a shaded zone. 



5.3.2 Asynchronous XNOR and XOR gate 

Exploiting some features of the interaction between gl glider and ol oscillator 



(Fig. 24) we can implement an asynchronous device which calculates XNOR 



and XOR operation at once (Fig. 251. Such a gate is designed by a scheme 



similar to that outlined in [B]. Oscillator in position shaded by gray in Fig. 24 
represents logical value True and absence of the oscillator — value False 
of logical operation XNOR, exclusive NOR operation. An auxiliary oscillator 
(generated when oscillator in shade position is annihilated) represents results of 
XOR operation, exclusive OR operation. 

Assuming that signals on inputs x and y can be generated at any time (no 
synchronization) except not at the same time, we obtain the following dynamics 
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Figure 25: Scheme of XNOR and XOR gate. 



of the device (Fig. 25): 
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XNOR 


XOR 








ol 








gl 





ol 


gl 








ol 


gl 


gl 


ol 






where ol and gl stay for oscillator and glider, means absence of the objects. 

6 Discussion 

Findings discussed in the paper are based on computational experiments with 
the Diffusion Rule CA and, particularly, exhaustive search of mobile and sta- 
tionary localizations emerged in spatio-temporal dynamics of the automaton. 




Figure 26: Simulating a luminescence pattern with the Diffusion Rule. 
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Amongst known 2D CA supporting localizations, the Diffusion Rule CA is 
the minimal model because cell-state transitions depend not on intervals of 'cell 
sensitivity' but on singletons, i.e. transition 0^1 occurs if there is exactly 
two neighbors in state 1 , and transition 1 ^ 1 if there exactly seven neighbors 
in state 1. Moreover, we are not aware of any other CA which exhibits so 
large variety of mobile localizations (gliders) and high diversity of outcomes 
of collisions between mobile and/or stationary localizations. Despite of trying 
to undertake exhaustive study of localization dynamics we nevertheless missed 
several important points that could form objectives of future studies. 

A simple but interesting physical simulation was made setting a reaction like 
luminescence. Using packages of both diagonal lines of 50 cells everyone. The 
luminescence phenomenon was obtained over the evolution in densities of cells in 



small intervals as illustrated the Fig. 26 The final state is dominated by blinkers 



or oscillators. This simulations are developed by the group of researchers iGEM- 
Mexico at the MIT0 

It is not a trivial problem to find large stable patterns constructed with big 
complex structures. Cell colonies damaged by a virus is one of such configura- 



tions (Fig. 27 1 



Diffusion Rule can simulate the evolution of a cell like the elemental cellular 
automata (ECA) Rule 18 or Rule 90. For example, we take a diagonal line with 
503 cells in the initial condition. During the evolution original line is multiplied 
in lines less and less small producing oscillators. Finally, the evolution space is 
dominated by oscillators that represent exactly a cell alive in ID case. However, 



there is generated a second evolution of the same type as illustrated in Fig. [28 
In this example, it took 512 steps of evolution to reach the configuration of 7,596 
living cells. 

The existence of stable configurations seems difficult to find in a rule which 
is generally chaotically producing super-nova explosions. Nevertheless, we have 
an example where four glider guns are synchronized to annihilate the gliders. In 
this case, two gl gliders and two g4 gliders come into quadruple collision shown 
in Fig. [29| 

We envisage that important open problems to be solved include implemen- 
tation of quasi-chemical reactions between gliders, studies of grammars derived 
and implementation of a full effective decision procedures based on glider col- 
lisions. It will be also worth to demonstrate intrinsic universality and self- 
reproduction. Another project would be to use de Bruijn diagrams [23] to check 
if there are any still undiscovered gliders with velocity one or still life con- 
figurations, and besides to use algorithms specialized in automatic search for 
complicated or big gliders [T5], oscillators, glider guns or more configurations. 
Also we are planning to make an exploration of the cluster of semi-totalistic 
rules originated by the Diffusion Rule[^ Finally, the last but not least open 
problem is to decide if all types of gliders can be constructed in collision with 
other gliders, a closure property with respect to set of gliders. 



http : //www. f enomec .unam.mx/pablo/ igem/ 

9 h-hml 

J 



^'http : //uncomp .uwe . ac .uk/genaro/Lif e_dc22 .html 
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Figure 27: Virus propagation in the Diffusion Rule and three mutations in 
same initial condition. Evolution rules and data are given below snapshots. 



26 



Figure 28: Simulating the evolution of a cell of the EGA Rule 18 or Rule 90 
with the Diffusion Rule. 
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Figure 29: Simultaneous annihilation across of four particles produced by four 
glider guns in the Diffusion Rule. The evolution shows a global configuration in 
36 steps with 118 live cells. 
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